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The percentage of substitutional doping of magnetic atoms (Mn) in group-IV-based dilute mag- 
netic semiconductors (DMS) can be increased by co-doping with another conventional electronic 
dopant (e-dopant), as demonstrated from first-principles calculations recently [Zhu et aL, Phys. 
Rev. Lett. 100, 027205 (2008)]. Here, we report extensive theoretical investigations of the kinetic 
and thermodynamic characteristics of several co-doped systems including bulk Si and Ge as hosts 
and various group-Ill and group- V e-dopants. The main findings are as follows: The n-p pairing 
of n-type e-dopants with p-type substitutional Mn is energetically stable in bulk Ge and Si. Mn 
atoms move from interstitial sites to substitutional sites easier (with lower kinetic barriers) in the 
presence of a neighboring n-type e-dopant. Magnetic coupling between two Mn atoms in bulk Ge 
oscillates between positive (ferromagnetic) and negative (antiferromagnetic) values with increasing 
Mn-Mn distance, but in Mn/As co-doped Ge the coupling parameter remains positive at all dis- 
tances beyond nearest-neighbors and this qualitative difference does not change with the doping 
level. For Mn doped Si, all coupling values except for the nearest neighbor one are positive and do 
not change much upon co-doping. We find an unconventional magnetic anisotropy in the co-doped 
system, that is, the dependence of magnetic coupling on the relative positions of the magnetic ions 
and their neighboring e-dopants. We map the calculated magnetic coupling to a classical Heisenberg 
model and employ Monte Carlo simulations to estimate the Curie temperature (Tc). We find that 
in Mn doped Ge no ferromagnetic order exists for Mn concentrations ranging from 3.13% to 6%. 
Instead, a spin-glass phase transition occurs at ^-^5 K at 5% Mn doping. For Mn/As co-doped Ge, 
Tc increases nearly linearly with the Mn concentration and reaches 264 K at 5% Mn doping. 



PACS numbers: 75.50.Pp, 66.30.Jt, 75.30.Hx 
I. INTRODUCTION 

Dilute magnetic semiconductors (DMS) have attracted 
much interest in the condensed matter community not 
only because of their promising application in spintronic 
devices^ ^, but also because of the many new and im- 
portant theoretical issues which arise from the study of 
this unique class of disordered magnetic system ^ l ^ l ^ l '^ l 
As for specific materials, besides the most extensively 
studied (III,Mn)V system^^, Mn doped group-IV semi- 
conductors such as Ge and Si also show promise for 
real applicationPISIIIIIlIinilllill^ ^^^^^ realize this 

promise, a Curie temperature comparable to room tem- 
perature or higher is required. Both theory and experi- 
ment indicate that the Curie temperature of the above- 
mentioned materials is exceptionally sensit ive to the r atio 
of interstitial to substitutional Mn atoms^^^E^HHni. in 
(III,Mn)V as well as (Mn,IV) systems, substitutional Mn 
atoms act as acceptors and provide holes which, accord- 
ing to current understanding, are the mediator of mag- 
netic interactions between magnetic moments in these 
materials. Interstitial Mn atomJ^ are identified to be 



donors, and tend to compensate the holes and magnetic 
moments induced by the substitutional MeF. Further- 
more, though annealing is an effective way to decrease 
the percentage of interstitial Mn while keeping the homo- 
geneity in (Ga,Mn)As^^^ it is less useful for Mua^Gei-a^ 
and Mn^Sii_^,'221^^1^^ which makes it very difficult to 
get high quality samples of these materials using conven- 
tional methods. 

In our recent work^^, a novel way to enhance the sub- 
stitutional doping of Mn in Ge and Si was proposed. 
In this method, an additional conventional electronic 
dopant (e-dopant) such as As or P is introduced in the 
doping process. Using first-principles electronic structure 
calculations, we were able to show that the co-doping ap- 
proach can substantially lower the energy of Mn atoms 
at substitutional sites relative to that at interstitial sites, 
as well as the energy barrier which the Mn atoms have to 
overcome in order to be incorporated into substitutional 
sites. In addition, the assisting e-dopant enhances the 
magnetic coupling between substitutional Mn atoms. A 
new type of magnetic anisotropy was also found, which 
depends on the proximity of the assisting e-dopant to a 
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Mn dopant, rather than the direction of magnetic mo- 
ments relative to the lattice direction of the host or to 
other moments. 

In this paper, we present a detailed ab initio investiga- 
tion of this novel approach by analyzing the kinetic and 
thermodynamic issues related to the stability of various 
dopant-host combinations. We then calculate the mag- 
netic coupling between two Mn atoms in bulk Ge and Si, 
and find that in Ge, the coupling oscillates between pos- 
itive (ferromagnetic) and negative (antiferromagnetic) 
values with the Mn-Mn distance. But in Mn/As co- 
doped Ge the coupling parameter remains positive at all 
distances beyond nearest-neighbors, and this qualitative 
difference does not change with the doping level. For 
Mn doped Si, all the couplings except for the nearest 
neighbor one are positive and do not change much upon 
co-doping. We also carry out Monte Carlo simulations 
to obtain the Curie temperatures of the co-doped mate- 
rials. We find that in Mn doped Ge no ferromagnetic 
order exists for Mn concentrations ranging from 3.13% 
to 6%. Instead, a spin-glass phase transition occurs at 
-5 K at 5% Mn doping. For Mn/As co-doped Ge, 
increases nearly linearly with the Mn concentration and 
reaches 264 K at 5% Mn doping. 

The paper is organized as follows: In section [ll] we 
present the computational methodology, including details 
of our ab initio treatment and the Monte Carlo simu- 
lations. The main ab initio results are given in Sec |III[ 
where the kinetic and energetic properties of various com- 
binations of host materials (Ge, Si) and assisting dopants 
(As, P, Al, Ga) are investigated. Magnetic interactions 
in As co-doped Mua^Gei-.^ are discussed in Sec|IV| In 



Sec|V|the ab initio results of magnetic coupling are used 
to find the transition temperature Tc of Mn/As co-doped 
Ge. The discussion and summary are provided in the last 
two sections. 



II. METHODS 



on Mn concentration. Specifically, we used a 3 x 3 x 3 
supercell, which corresponds to 216 atoms, and with one 
of them replaced by a Mn the concentration is 0.46%. In 
each calculation of the magnetic coupling between Mn 
atoms, two Mn atoms are placed in a 3 x 3 x 3 super- 
cell, corresponding to a 0.926% Mn concentration. We 
also selectively use a 2 x 2 x 2 supercell for the mag- 
netic coupling with two Mn atoms in the supercell, cor- 
responding to 3.125% Mn, for comparison. This setup is 
similar to prev ious studie s of Mn-Mn interactions in pure 
semiconductor j^ ^'^ ' ^^ ' ^^ i 

A uniform 4x4x4 (2x2x2) mesh, including the 
F point (0, 0, 0), is chosen for Brillouin Zone sam- 
pling in the 2x2x2 (3x3x3) supercell. Optimized 
atomic geometries are obtained when the forces on all 
the unconstrained atoms are smaller in magnitude than 
0.01 eV/A. The "climbing image Nudged Elastic Band" 
(NEB) methocF^is used to locate the transition state ge- 
ometries for the calculation of activation energy barriers. 
Typically four slab replicas between the initial and final 
geometries are enough to produce a smooth minimum 
energy path. 

For the Monte Carlo simulations we use the Metropo- 
lis algorithm^^ and the magnetic energy of the system is 
calculated using the classical Heisenberg model, in which 
each magnetic ion is treated as a classical moment and 
is placed at a randomly chosen site of the supercell. The 
magnetic coupling parameters are extracted from ab ini- 
tio results of the energy difference between parallel and 
antiparallel spin configurations of two Mn moments at 
different separations. At each temperature we use 50000 
Monte Carlo steps per moment for the system to relax, 
and calculate the thermal average in the following 50000 
steps. To determine the Curie temperature, we adopt the 
fourth order cumulant crossing metho d bas ed on finite- 
size scaling theory proposed by BindeiE^EH. In applying 
this method we choose three supercell sizes: 8x8x8, 
10 X 10 X 10, 12 X 12 X 12, and 40 configurations in each 
case for averaging. 



Our spin-polarized first-principles calculations are car- 
ried out using the Vienna ab initio simulation package 
(VASP)^'', a density functional theory appro ach u sing 
the projector augmented wave (PAW) methocP^ and 
the generalized gradient approximation (PBE-GGA)I^ 
for exchange-correlation. A default plane- wave energy 
cutoff of 269.9 eV is consistently used in all Mn calcula- 
tions. These choices produce a bulk Ge and Si lattice con- 
stants of 5.78 A (experimental value^^ 5.66 A) and 5.47 
A respectively (experimental value^^ 5.43 A). 

In our calculations of the co-doping process the su- 
percell size is chosen to be a 2 x 2 x 2 multiple of the 
conventional cubic cell of the diamond lattice which con- 
tains 8 atoms. Hence, there are 64 atoms in one supercell, 
and one of them is replaced by an Mn atom, correspond- 
ing to 1.56% Mn concentration, comparable to what was 
achieved experiment all}/ ^ 34|35 i^ Different supercell sizes 
were used to study the dependence of calculated results 



III. AB INITIO STUDY OF THE CO-DOPING 
PROCESSES 

A. Study on intrinsic (Mn,IV) without additional 
dopants 

We first consider the equilibrium structure of a single 
Mn dopant atom in bulk Si and Ge, and address the diffi- 
culty of lowering the percentage of interstitial Mn impu- 
rities. A complete understanding of the microscopic dop- 
ing process requires detailed knowledge of the energetics 
as well as the kinetics of dopants in the host materiaP^. 
In fact, an in-depth understanding of the growth kinet- 
ics is particularly important, because DMS systems are 
typically in a metastable state, since they are grown by 
co-doping the magnetic dopants and host semiconductor 
atoms using molecular beam epitaxy under nonequilib- 
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FIG. 1: (color online) Different Mn sites in bulk Si or Ge: 
(a) Mn at a substitutional site, (b) Mn at an interstitial site, 
(c) Final state of an interstitial Mn kicking out a neighboring 
host atom to an interstitial site and occupying the left-behind 
substitutional one. 

rium conditiQn j^ l ^'^ 1 ^^ 1 ^^ 4 

To address these issues, we calculate the relative for- 
mation energy of a substitutional Mn(Fig. [TJa)) and in- 
terstitial Mn(Fig. [ijb)) atom in Ge and Si separately, 
which is defined as: 

AEi = {Esubst + 

where iihost is the host material's chemical potential. For 
Ge, our calculation gives A£^i = — 0.63eV. Thus, in Ge 
the substitutional sites have a relatively lower energy and 
are preferred by Mn atoms. However, for Si the opposite 
is true and A^i = +0.58. This reversed site preferencd^ 
makes it extremely hard to achieve experimentally even 
a nominal concentration of substitutional Mn in silicon. 

We next consider kinetic aspects of the Mn doping pro- 
cess. In order to get a high ratio of substitutional to in- 
terstitial Mn, the process that an interstitial Mn kicks 
out a host atom and becomes substitutional must take 
place more often than the reverse process. Accordingly, 
we calculate the energy difference between the initial (in- 
terstitial, ([T])(b)) and final (substitutional, ([l])(c)) states 
of this process: /S.E = Ef — Ei, and the energy barriers 
Ea and e'^ for the reverse process. Our calculation shows 
that for both Ge and Si AE is positive (0.82 eV and 2.03 
eV respectively). This energy cost for the transition from 
initial to final state defines the lower bound of the activa- 
tion energy barrier for the exchange process, which must 
be lower than ~ 0.8 eV for efficient incorporation ( with a 
standard attempt frequency 10^^ sec~^). Moreover, the 
actual energy barrier Sa in either case (1.12eV for Ge 
and >2eV for Si) is higher than the barrier of the reverse 
process, which is calculated as e'^=ea-AE^ with the latter 
being lower than 0.8 eV, further facilitating the reverse 
processes. Thus, kinetically Mn is more stable at inter- 
stitial sites rather than at substitutional sites in both Ge 
or Si. In the following section we will address the issue 
of doping Mn together with another n- (P, As) or p-type 



(Al, Ga) conventional e-dopant in order to explore how 
the assisting dopants influence this site preference both 
energetically and kinetically. 



B. Energetic and kinetic study on the co-doped 
systems 

Substitutional Mn in Ge is a p-type double acceptoi^^. 
Our proposal for a co-doping mechanism is based on the 
fact that the electrostatic interaction between a n-type 
and a p-type dopant in a semiconductor is attractive be- 
cause of their different charge states (see below). Thus 
an n-type e-dopant may help to stabilize substitutional 
Mn atoms. 

We start by noting that in Ge or Si there are two 
kinds of interstitial sites: the hexagonal interstitial site 
Ih^ which has six nearest neighbors, and the tetrahedral 
site It with four nearest neighbors. Using first-principles 
calculations we find that in n-type doped Ge and Si, the 
energy of a Mn sitting at the Ih site is different from 
that at the It site. For P, As and Sb doped Ge, the en- 
ergy differences are 0.14, 0.09 and 0.04eV, respectively, 
where a positive sign means the Ih occupation has lower 
energy and is preferred. In the case of either n- or p- 
doped Si as well as p-doped Ge, It is preferred to Ih- 
This dopant dependent preference can be qualitatively 
explained by the local strain effect. Namely, a Mn atom 
and an n-type dopant favor a relatively short bonding 
distance, which is accommodated by Mn occupying the 
Ih site rather than the It site in Ge (the Ih site has 
a shorter distance to its nearest neighbors than the It 
site). To show that this is indeed the case, we reduce 
the lattice constant of Ge to the value of Si and calculate 
the energy difference again. Then the results show that 
the preference for Mn is changed to the It site, because 
in this case the distance between the Ih Mn and n-type 
dopant becomes too short (compressive), whereas at the 
It site the Mn/ n-type dopant bond length is close to its 
optimal value. We have also checked to confirm that if 
we increase the lattice constant of Si to that of Ge, the 
preference for Mn is changed to the Ih site for the n-type 
doped systems. 

In the following we examine two possible kinetic pro- 
cesses of an interstitial Mn atom becoming substitutional. 
These processes share the same initial state with Mn oc- 
cupying either the Ih or It sites with a neighboring n- 
type or p-type e-dopant. From our calculation of the to- 
tal energy of a Ge-supercell with an interstitial Mn and 
a substitutional e-dopant as a function of their separa- 
tion, shown in Fig. |2j we find that shorter separation is 
energetically preferred. Thus, the choice of neighboring 
Mn/e-dopant pair configuration is reasonable. 

In the first process, denoted as Process I, we consider 
an interstitial Mn directly exchanging position with its 
substitutional e-dopant neighbor. In the final state, the 
e-dopant is pushed to an adjacent interstitial site and the 
Mn atom moves into the substitutional site left behind. 
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FIG. 2: Calculated relative total energy as a function of 
the distance between an interstitial Mn and a substitutional 
dopant in bulk Ge. 



FIG. 4: Calculated relative total energy as a function of the 
distance between a substitutional Mn and an interstitial n- 
type dopant in bulk Ge. 



as shown in Fig. [3] Table |l] summarizes the calculated en- 
ergy differences between the final and initial states 
for n-type and p-type dopants in Si and Ge. We find that 
only the P or As doped Ge (with A£:=0.33 eV and 0.42 
eV respectively) can fulfill the requirement that l\E <0.8 
eV. However, further examination of the activation en- 
ergy barriers for incorporation in these two cases gives 
£0= 0.88 eV and 0.98 eV respectively, which means this 
process is unlikely to happen in both cases. Moreover, 
the reverse processes with e'^-= 0.55 eV and 0.56 eV, re- 
spectively, are more likely to occur. 

Nevertheless, there is one possibility for the Mn atom 
to stay at the substitutional site, that is, the kicked-out 
dopant atom diffuses away rapidly so that the reverse 
process cannot happen. This is ruled out by our calcula- 
tion of the energy of a Mn/e-dopant pair as a function of 
their separation, shown in Fig|4j which shows that the e- 
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FIG. 3: (color online) Atomic structures and schematic en- 
ergy profiles for Process I: (a) Initial state with Mn in the 
It position (except for n-type doped Ge); (b) Final state; (c) 
Initial state for n-type doped Ge with Mn at the Ih position. 



dopant cannot diffuse away because the energy increases 
with increasing separation. 

We then consider a different process (Process II) , which 
starts from the same initial configuration as in Process 
I, but instead of exchanging with the e-dopant, the Mn 
atom now pushes out a host atom next to the e-dopant to 
an interstitial site, and then occupies the substitutional 
site left behind. The final state is shown in Fig. |5j in 
which the kicked-out interstitial host atom, the substitu- 
tional Mn and the e-dopant are nearly collinear. The cal- 
culated and Sa for various n-type and p-type dopants 
in Si and Ge is also shown in Table [l[ For n-type doped 
Ge these values are substantially lower than in Process I, 
and considerably below the threshold value of 0.8 eV. For 
P and As doped Ge, is actually quite low. Further- 
more, the activation barriers Sa for all the three n-type 
dopants is lees than 0.4eV. Qualitatively, this substantial 
change in the energetic and kinetic characters originates 
from the electrostatic attraction between the Mn atom, 
which behaves like a p-type dopant, and n-type dopants. 
Therefore, Process II, leading to substitutional Mn atoms 
proximate to n-type dopants, is more likely to happen in 
reality. 

One issue that arises at this stage is whether the fi- 
nal state is thermodynamically stable. To address this 
question, we calculate the energy difference between in- 
terstitial Mn and substitutional Mn defined as 

A£^2 = {Epair + l^host) ~ Winter - (2) 

Here Einter is the total energy of a Mn/e-dopant pair, 
with the Mn sitting at an interstitial site, while Epair 
is that with the Mn occupying a substitutional site. 
The calculated interstitial-substitutional energy differ- 
ence AE2 is shown in Table [llj Compared to the results 
without n-tyep dopants in Sec. |III A[ the substitutional 
Mn in Ge becomes much more stable with the neighbor- 
ing n-type dopant. Moreover, the site preference of Mn 
in Si is reversed from interstitial to substitutional. 
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TABLE I: Calculated energy differences — Ef — Ei (in eV) between the final and initial states of Process I and Process 
II, illustrated in Fig.l. Sa (in eV) is the activation energy for a transition from the initial to final state. Results highlighted 
in bold correspond to processes for which or Sa or both are < 0.8 eV. All the results are for Mn concentration of 1.56%; 
results for selected cases with Mn concentration of 0.46% are given in brackets. 
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FIG. 5: (color online) Atomic structures and schematic en- 
ergy profiles of Process II: (a) Initial state with Mn in the 
It position (except for n-type doped Ge); (b) Final state; (c) 
Initial state for n-type doped Ge, with Mn at the It position. 



We next calculate the total energy of a Ge-supercell 
doped by a substitutional Mn/e-dopant pair at different 
separations. The trend of the total energy with increas- 
ing distance between the Mn atom and the e-dopant is 
shown in Fig. |6] The interaction between Mn and e- 
dopant is attractive for n-type e-dopants (P and As) and 
repulsive for p-type e-dopants (Al and Ga). This suggests 
that the picture of electrostatic interaction between Mn 
and e-dopants that we proposed at the beginning of this 



TABLE II: Relative formation energy of substitutional and 
interstitial Mn in the presence of a neighboring substitutional 
n-type dopant, defined as: A^2 = {E pair + IJihost) - E inter (in 
eV). Negative values indicate higher stability of the substitu- 
tional configuration over the interstitial. The relative energy 
of substitutional and interstitial Mn in pure Si or Ge are in- 
cluded for comparison. 
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section is valid. 

Finally we note that the energy differences between 
the initial and final states depend on the Mn concentra- 
tion, as illustrated in Table I. The calculated AE values 
at the 0.46% Mn concentration are larger than those at 
the 1.56% concentration, but for the important cases of 
n-type dopants in Ge, these energy differences are still 
much lower than the threshold of ~ 0.8 eV. This rela- 
tively strong dependence is not due to constant volume 
calculations, because it is also observed when the super- 
cell volume is fully relaxed. Instead, it is caused by the 
interaction between the Mn atoms in adjacent supercells. 
We stress that the qualitative picture that the n-type 
dopants facilitate substitutional incorporation of Mn is 
valid for all the experimentally accessible Mn concentra- 
tions considered here. 

In short, we have shown that in the presence of a neigh- 
boring n-type dopant, the substitutional sites are ener- 
getically preferred by Mn atoms to interstitial sites and 
are kinetically accessible. In the following sections we 
will turn to study the electronic and magnetic properties 
of this new n-p co-doped system. 
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FIG. 6: Relative total energy as a function of the distance 
between a substitutional Mn and a substitutional e-dopant in 
bulk Ge. 
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IV. ELECTRONIC AND MAGNETIC 
PROPERTIES OF THE CO-DOPED SYSTEMS 

A. Electronic structure 

The electronic properties of the Mn/e-dopant co-doped 
systems are conveniently presented through the calcu- 
lated density of states (DOS). Fig. [t] and Fig. [s] show the 
total DOS and local DOS for the substitutional Mn and 
e-dopants in Ge and Si, respectively. Several important 
features emerge: 

(1) Mn doped Ge or Si are all half metals, regardless of 
the existence of e-dopants like As or P, which means the 
value of the total magnetic moment per Mn atom is in- 
teger. 

(2) From the figures it can be determined that the mo- 
ment per Mn is Sjhb in pure Ge or Si, and after 
co-doping with another e-dopant. The importance of this 
finding, namely, co-doping can actually increase the mag- 
netic moment of Mn, will be discussed in the next sub- 



section. 

(3) The local DOS for Mn is broadened to the whole 
range of the host valence band, indicating that there is 
strong hybridization between the Mn d-state and the va- 
lence p-state of the host semiconductor. 

(4) The local DOS of the additional e-dopant is negli- 
gible, meaning that the states it contributes are mostly 
delocalized, so that its most evident infiuence on the to- 
tal DOS is simply to shift the Fermi energy to a higher 
value, as is expected for a regular non-magnetic dopant. 



B. Magnetic properties 

At first sight, the additional n-type dopant may nega- 
tively infiuence the strength of the magnetic interaction 
betwee n Mn atoms, because of the compensation of hole 
carrierJ^J^. This argument may not be true for the fol- 
lowing reason: We can write the magnetic interaction 
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FIG. 7: (color online) The spin-resolved DOS of (a) Mn 
doped Ge, and (b) a Mn/As pair doped Ge. Projected DOS 
of Mn 3d states and As are given as red and blue dashed lines, 
respectively. The DOS for bulk Ge is shown for comparison. 



FIG. 8: (color online) The spin-resolved DOS of (a) Mn 
doped Si, and (b) a Mn/P pair doped Si. Projected DOS of 
Mn 3d states and P are given as red and blue dashed lines, 
respectively. The DOS for bulk Si is shown for comparison. 
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energy between two Mn ions as: 



n\—Mn2 



(3) 



where Jeff denotes the effective magnetic couphng 
strength and Smui^ SMn2 represent the local magnetic 
moments associated with the Mn atoms; even if Jg// were 
to decrease because of the carrier compensation effect, 
since the influence of the e-dopant on the local moments 
of Mn is positive as mentioned in point (2) of the previ- 
ous subsection, it is still possible that the enhancement 
of Mn moments by the e-dopants outweighs its negative 
influence on Jeff- 

To check whether this is the case, we first resort to 
direct ah initio calculation of the magnetic coupling en- 
ergy of a Mn-Mn pair with different separations, which 
can be represented by the total energy difference 
between the antiferromagneti c (AFM) a nd the ferromag- 
netic (FM) states of the pai j^^ ' ^^'^^ ' ^ . In the present 
case, each Mn atom has an n-type dopant neighbor, 
which leads to more spatial configurations with the same 
Mn-Mn distances. To be precise, a substitutional Mn 
atom has four nearest neighbors, that is, four possible 
sites for the n-type dopant atom, and thus there are 16 
possible configurations for a given Mn-Mn distance. The 
number of nonequivalent configurations for each of 12 
Mn-Mn separations in the range 2.4 - 9.5 A in Si (2.5 - 
10.0 A in Ge) is 2, 6, 16, 4, 7, 16, 4, 7, 16, 10, 7 and 5 
respectively, with increasing distance. 

With these considerations, our results for Mn/As co- 
doped Ge and Mn/P co-doped Si are shown in Fig. |9] 
We first note that in the case of two Mn atoms in pure 
Ge ((a) and (b) of Fig. |9|, the behavior of the AFM- 
FM energy difference is oscillatory between posi- 
tive and negative values as a function of distance (but 
monotonically decreasing along different directions, see 
below). In contrast, the average interaction between the 
two Mn/As pairs in Ge always favors FM coupling except 
at the nearest-neighbor Mn-Mn distance, and this char- 
acteristics does not change with doping level (compare 
Fig.[9];a) of 3.125% Mn and Fig.[9];b) of 0.926% Mn). In 
the case of Si as host, Mn atoms favor FM coupling ex- 
cept for nearest neighbor distance, and this feature does 
not change upon co-doping. 

Though usually it is assumed that the hole-mediated 
magnetic interaction in dilute magnetic semiconductors 
is RKKY-like, we find that the FM-AFM oscillation dis- 
played in Fig. [9] should not be treated as a manifesta- 



tion of the RKKY interactioiPE^^^. Fig. IToFa) shows 



the magnetic coupling between Mn ions along different 
lattice directions, which is similar to the design in Ma- 
hadevan's worlP^, where it is evident that the oscilla- 
tory behavior is replaced by monotonic decrease in mag- 
nitude. On the other hand, for the doping levels con- 
sidered here, the period of RKKY oscillation is much 
larger than the lattice const ant^^'^^, as in the case of 
GaAs. Thus, the oscillation here is merely due to mag- 
netic anisotropy, rather than a manifestation of RKKY- 
type interaction. In Fig.ITOfb), we plot the coupling along 
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FIG. 9: (color online) Total energy difference between AFM 
and FM states of two Mn ions versus Mn-Mn separation for 
two Mn/As pairs doped Ge at (a) 3.125% and (b) 0.939% Mn 
concentration; (c) for two Mn/P pairs doped Si, represented 
by small orange dots. The large red dots are averages over 
the small orange dots for a given Mn-Mn distance. For com- 
parison, the results for the systems doped with only two Mn 
impurities are shown as blue diamonds. 
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FIG. 10: (color online) Magnetic coupling along different 
crystal directions. Black squares are for the (111) direction, 
red dots for the (110) direction and blue triangles for the (100) 
direction. Magnetic coupling along a close-packed chain is 
illustrated in the inset. 



a close-packed atom chain, which also shows monotonic 
decrease with distance except for the nearest neighbor 
value. This finding leads us to suggest that Mn ions in 
Ge are magnetically coupled through some paths consist- 
ing of covalently-bonded Ge atoms, a hypothesis which 
deserves to be checked by more detailed investigations. 

There are two other important issues revealed in Fig.|9] 
First, the dispersive values of AE at a given distance 
in the co-doped case show a new kind of magnetic 
anisotropy. This anisotropy is conceptually different from 
the magnetic anisotropy that is typically discussed in the 
literature. The traditional definition refers to the cases in 
which the magnetization of a system exhibits anisotropy 
when the magnetic moment is polarized along different 
crystalline directions'^, or when the coupling between 
two magnetic dopants is anisotropic depending on their 
relative orientation in the host material. In contrast, 
here the two magnetic dopants are fixed in space, and the 
magnetic anisotropy is caused by the relative positions of 
the two n-type e-dopants surrounding the magnetic im- 
purities. Second, from the presence of the n-type dopant. 



the FM interaction between two magnetic atoms on the 
whole preserves its magnitude rather than being substan- 
tially weakened. Thus, the influence of e-dopants on the 
magnetic properties of the whole system is not simply 
a weakening of the magnetic coupling by decreasing the 
number of interaction mediators. 

Deeper understanding of the above observations re- 
quires a careful examination of the microscopic coupling 
mechanism. To this end, we consider three representa- 
tive configurations, all with the same Mn-Mn separation 
fixed, equal to the next nearest neighbor distance in the 
Ge matrix: (a), a Mn-Mn pair in pure Ge as the reference 
structure; (b) and (c), a Mn/e-dopant-Mn/e-dopant pair 
with the strongest and weakest magnetic couplings, re- 
spectively. For the reference case shown in Fig. llja), the 
two Mn atoms share a Ge atom as their nearest neighbor. 
When the Mn pair is ferromagnetically coupled, the spin 
density in the plane containing the two Mn atoms and 
their mutual Ge neighbor (the (110) plane indicated 
in Fig. [nja)) is plotted in Fig. [iTJb). The red (blue) 
area represents spin up (down) density. The large local 
magnetic moments of Mn induce spin polarization on the 
nearby non-magnetic Ge atoms, which are antiferromag- 
netically coupled with the Mn atoms. 

The corresponding plots for case (b) are shown in 
Fig. [TTJc) and Fig.jlTJd). In this case, the two As atoms 
are not in the (1 1 0) plane, and the two Mn atoms still 
have the same Ge atom as their mutual nearest neighbor. 
Furthermore, the local magnetic moment of the bridging 
Ge atom shows little change, indicating that Jeff essen- 
tially stays the same. To show the effect of As doping, 
we plot in Fig. 11 ^f) the spin density on the plane con- 
taining the two Mn and one As atom (the (112) plane in 
Fig. [TTJe)). Here, As acts as a donor helping to compen- 
sate the holes introduced by its neighboring Mn, resulting 
in an increased local magnetic moment on each Mn atom, 
Smu (3.60/iB^ 4.00/iB). Therefore, the overall magnetic 
coupling between the two Mn atoms is enhanced relative 
to the pure Ge case. 

For case (c), the corresponding plots are shown in 
Fig. [njg) and Fig. [iTJh). In this case, the two As and 
two Mn atoms are both in the (110) plane, with one As 
replacing the mutual nearest Ge neighbor of the two Mn 
atoms. Similar to case (b), here both 5'Mniand6'Mn2 are 
also increased (5'Mni'3.60/iB^ 3.87/iB; 5'Mn2-3.60/iB^ 
4.02/i^; the asymmetry in the increase is caused by 
the asymmetric locations of the two As atoms). How- 
ever, because the local magnetic moment of the bridg- 
ing atom is substantially decreased from that of case (a) 
(Ge:— 0.16/iB-^ As:— O.OS/i^), the corresponding Jg// is 
also significantly weakened, leading to an overall weak- 
ened magnetic coupling between the two Mn atoms rela- 
tive to the pure Ge case. 

Summarizing, As as n-type dopant can enhance the 
local magnetic moments of neighboring Mn atoms, but 
itself is weakly spin polarized (much weaker than Ge). 
Therefore, if As serves as the bridging atom between two 
Mn atoms, the global magnetic coupling will be weak- 
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FIG. 11: (color online) The atomic structures and spin den- 
sity plots of three representative configurations of two Mn 
TOMS, (a) and (b), and two Mn/As pairs, (c)-(h), doped 
Ge. In all structures, the two Mn atoms are fixed at the next 
nearest neighbor distance. The spin density plots are taken 
on the green plane as depicted in the structures on the left. 
The red and blue contours represent the two different spin 
components, (c) and (e) correspond to the configuration with 
the strongest magnetic coupling between two Mn ions, and 
(g) the weakest magnetic coupling. 



ened. If As is located so as to only enhance the magnetic 
moment of Mn, with a Ge atom still bridging the Mn-Mn 
coupling, then the global magnetic coupling will be en- 
hanced. This conclusion is further confirmed by checking 
other Mn-Mn distances. 



V. CURIE TEMPERATURE 

To study the macroscopic magnetic properties of the 
co-doped DMS materials using our ab initio results, we 
turn to the classical Heisenberg model: 



J' 



(4) 



where Jij is the magnetic coupling constant between mo- 
ment i and j, and Gi is a unit vector representing the 
direction of spin i. Then the AFM-FM energy difference 
calculated in previous section is given by: 



l\E — Eafm — EpM = 4Ji: 



(5) 



with 1 and 2 the indices of the two moments in the super- 
cell. With given coupling parameters, we then use Monte 
Carlo simulations to address the statistical mechanics of 
the DMS systems at finite temperatures. To eliminate fi- 
nite size effects, the cumulant crossing method^^ is used 
to determine the Curie temperature. This two-step ap- 
proach has the distinct advantage over the ordinary mean 
field approach, that both disorder and percola tion e ffects 
are naturally and precisely taken into accounlP^^. 

A subtle issue in the present case is the following: 
Since real interactions between magnetic atoms in DMS 
have a built-in multiatom nature, in an optimal Heisen- 
berg description the coupling parameters must depend 
on the system geometry. This is difficult, if not entirely 
impossibl d^"*^ * ^^^ to address, because of the very large num- 
ber of possible configurations in a macroscopic system 
and some approximations are necessar y. 

The supercell ab initio approactP^^^^^^I em- 
ployed here assumes only one approximation, the pair- 
superposition approximation, which means that the in- 
teraction is exclusively pairwise and can be added inde- 
pendently to get the total interaction. Though this may 
not hold at high concentrations of magnetic moments, we 
claim that it should be a reasonable approximation at the 
low concentrations we considered (3.13%-6%), where the 
average distance between two Mn atoms, calculated by 



d- 



327rx 



(6) 



with X the concentration and a the lattice constant, 
ranges from 1.97a to 1.58a, or 11.39A to 9.14A in the 
case of Mn doped Ge. Considering the bond length di^ond 
in Ge is only about 2.5A, a separation ~ 4:dbond is large 
enough for the system to be treated in this approxima- 
tion. 

Using the ab initio coupling parameters for Mua^Gei-a^, 
we first find that MC does not yield identifiable Tc up to 
X = 6% (see below). Nevertheless, after co-doping with 
As, MC shows that the system has high Tc, as summa- 
rized in Fig. [T2j in which we also include the results from 
the mean field approximation (MFA) using the formulsP^ 



1 2x 



(7) 



These results show that the MFA greatly ov erestim ates 
the Curie temperature, as established beford^^^. At 
X = 5%, Tc is evaluated to be 264K through MC, which 
is much higher than the 118K of 5% Mn doped GaAM 
At the 6% Mn concentration, MC gives a Tc higher than 
room temperature. Considering that x=6% is already a 
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FIG. 12: (color online) Comparison between Curie temper- 
atures calculated by the Monte Carlo approach (MC, black 
squares) and those obtained by the mean field approximation 
(MFA, red dots). 



relatively high concentration, we expect that the pair- 
superposition approximation may not be valid in this 
case. Arsenic doping can still be expected to dramati- 
cally change the magnetic properties of Mn doped Ge, 
namely, from no finite Tc to a potentially high Tc DMS 
material. 

The dependence of the Gurie temperature on Mn 
concentration, as obtained from the MG results, is al- 
most linear. This behavior is partly due to the pair- 
superposition approximation we used, meaning that the 
strength of magnetic coupling does not depend on Mn 
concentration. The only influence of concentration on 
Tc is the average number of magnetic impurity atoms on 
each coordination shell. Thus, after the configurational 
average, we expect that the dependence of Tc on x re- 
sembles the linear one obtained within the MFA (Eqn.[7|). 
Another reason for this linearity is that the concentra- 
tions we studied are higher than the magnetic percolation 
threshold of this systenP^. 

The presence of AFM couplings and the absence of Tc 
in the case of pure Mn doped Ge suggests the possibil- 
ity of a spin-glass ground state for this system. Jaeger 
et al^^ claimed that at low temperatures Mua^Gei-a^ ex- 
hibits spin-glass-like behavior and the critical tempera- 
ture of the spin-glass phase transition is 12K and 15K, 
for Mn concentrations x = 0.04 and x = 0.2, respectively. 
To examine whether this is the case, we first study the 
spin-spin correlation function of 5% Mn doped Ge, at 
T = O.Oli^. The result is shown in Fig. [Tsj along with a 
plot for Mn/As co-doped Ge, for comparison. The cor- 
relation function of Mn^^Gei-a; decays very fast with in- 
creasing distance and approaches to zero, indicating the 
absence of FM order even at low temperatures. Using 



FIG. 13: (color online) Spin-spin correlation function ob- 
tained from the Monte Carlo simulations. 



the spin-glass order parameter, defined as 

i 

and a similar cumulant crossing method^^ yields a tran- 
sition temperature ~5K, a value in semi-quantitative 
agreement with the results of Jaeger et a/.^. 



VI. DISCUSSION 

Mn^^Gei-a; has been attractive within the DMS com- 
munity because of its easy incorporation into the cur- 
rent semiconductor industry. The mechanism of va- 
lence hole mediated ferromagnetism for (Ga,Mn)As was 
proposed years ago ^ ' ^^ ' ^-^ ' and has been extensively ac- 
cepted ever since, but there is still no definitive theory for 
Mn^^Gei-J^. One reason, which is also one of the main 
points of this paper, is the difficulty of decreasing the per- 
centage of interstitial Mn dopants. The other important 
point is the hard-to-control inhomogeneity of this system, 
which has been realized only in recent years. The high 
Gurie temperature formerly reported in Mn^^Gei-a^^^^ 
is now thought to be due to the formation of Mn- 
rich regions in the host semiconductoi^ ^ \ For 

example, Mn-rich nanodots^^ and nanocolumn j^^ * ^^ -^ 
in Mn^^Gei-a; have been reported by many experi- 
mental groups, and later reproduced in Monte Carlo 
simulations^. 

Despite the seemingly unavoidable precipitation or 
spinodal decomposition^^ present during the growth of 
MiixGei-x samples, the study on homogeneously d oped 
Mn^^Gei-a; has never stopped. Work by Li et a/.EM] 
indicates that the long-range FM order in Mn^^Gei-a^ 
only exists at low temperatures (< 12K). Jaeger et alW^ 
claimed that even at low temperatures Mn^^Gei-a; shows 
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spin-glass-like behavior, and proposed that this is due 
to the intercluster frustration between FM Mn-rich clus- 
ters. Recently, Zeng et al}^ ^ using a newly developed 
subsurfactant epitaxy method, successfully grew cluster- 
free Mna^Gei_a^ samples with a Mn doping level of 0.25%. 
Surprisingly, this low doping level (by normal DMS stan- 
dards, where 1% to 5% is typical) led to a Curie temper- 
ature as high as over 400 K. 

The results in the present work provide a viewpoint 
that may resolve the seemingly conflicting experimental 
results discussed above. Specifically, we showed that the 
magnetic coupling between Mn ions in Mn^^Gei-a^ os- 
cillates between FM and AFM with increasing Mn-Mn 
distance and that homogeneous Mn^^Gei-a^ exhibits spin 
glass behavior. Thus, the FM order observed in experi- 
ments could be due to spatially ordered structures, which 
are formed due to precipitation or spinodal decomposi- 
tion. The high transition temperatures are expected be- 
cause of the large concentration of magnetic moments 
within the clusters. On the other hand, the AFM frus- 
tration in this case only manifest itself in the inter-cluster 
interaction, and thus leads to the spin-glass behavior at 
low temperatures. 

The unexpected high Tc in Zeng's work requires more 
discussion. Upon co-doping with n-type dopant As, the 
AFM coupling between Mn ions is absent, and a high 
Curie temperature emerges. We thus speculate that 
the high transition temperature in this case originates 
from this co-doping effect, and the unexpected n-type 
e-dopant here is most probably oxygen. Indeed, a re- 
cent study on the role of oxygen defects in Mn^^Gei-a; by 
Continenza and Profeta^^ supports this scenario, namely 
that oxygen acts as an n-type dopant and facilitates the 
substitutional Mn doping. It is also reasonable to expect 
a positive infiuence of oxygen on the Mn-Mn magnetic 
coupling, which, together with the possible existence of 
Mn-rich regions, can lead to a high Curie temperature. 

Recently, the works of KurodsP^ and Bonanni^^ 
demonstrated experimentally that the aggregation of 
magnetic ions in DMS systems can be controlled by mod- 
ifying the charge states of the magnetic dopants. This is 
in agreement with the spirit of our work, that is, charge 
states of impurity dopants play an important role in the 
growth kinetics of DMS materials and can lead to differ- 
ent structures with their own specific properties. 

Finally, this work suggests that the enhancement of 
substitutional Mn concentration in group-IV DMS can 
be achieved in epitaxial growth by co-depositing with 
the e-dopants. More specifically, this co-doping method 
can be integrated in the recently developed subsurfactant 
epitaxial growthP^, where pure Ge layers epitaxially grow 
on a Ge(lOO) substrate pre-covered with a submonolayer 
of Mn. During the growth process, the Mn atoms tend 
to diffuse upward to the subsurface layer, as predicted 
in a previous theoretical stud}' ''^. When the growth is 
slow enough, a small fraction of the Mn atoms can be 
trapped in substitutional sites, which leads to homoge- 
neous substitutional Mn doping with all the interstitial 



Mn fioating at the subsurface layer. However, the result- 
ing Mn concentration is still pretty low (0.25%). Here we 
propose that by co-depositing Ge with another n-type e- 
dopant, with very low depositing rates, the growth front 
could mimic the subsurfactant growth mode, but with 
more efficient substitutional trapping of Mn. The Mn 
trapping rate can be controlled by changing the concen- 
tration of the e-dopant. Experimental confirmation of 
this co-doping scheme is highly desirable. 



VII. CONCLUSION 

In conclusion, our ah initio DFT calculations show that 
in DMS materials additional n-type electronic dopants 
can serve to enhance the substitutional doping of p-type 
magnetic dopants such as Mn in the host group IV semi- 
conductors Si and Ge. The additional dopants suppress 
to a large extent the charge and magnetic-moment com- 
pensating effects from interstitial Mn, which is detrimen- 
tal to FM order. We calculate the magnetic coupling be- 
tween moments associated with Mn atoms using the en- 
ergy difference between parallel and antiparallel aligned 
pairs of Mn moments. We examined the unconventional 
magnetic anisotropy in Mn/As co-doped Ge, namely, the 
dependence of magnetic coupling on the relative posi- 
tions of magnetic ions and their neighboring assistant 
dopants. We find that the coupling oscillates between 
ferromagnetic (FM) and antiferromagnetic (AFM) with 
increasnig Mn-Mn distance in the Mn-doped Ge, whereas 
in As/Mn n-p co-doped Ge the coupling values at Mn- 
Mn separations up to the 12th coordination shell are all 
FM, except for the nearest-neighbor one. Moreover, we 
find that the FM-AFM oscillatory behavior in Mn^^Gei-a^ 
is due to anisotropy rather than being the result of a 
RKKY-type interaction. Our Monte Carlo simulations, 
using magnetic coupling parameters obtained from the 
ah initio calculations, indicate a high Curie temperature 
in Mn/As-Ge of 264K at 5% Mn doping. On the other 
hand, no FM order is observed in Mn^^Gei-a; (without 
co-doping) as Mn concentration ranges from 3.13% to 
6%. Thus, the homogeneously doped Mn^^Gei-a; is most 
likely a generic spin glass, with a spin-glass transition 
temperature of 5K at 5% doping, also obtained from our 
Monte Carlo simulations. Accordingly, we suggest that 
the high Curie temperature observed experimentally in 
MiixGei-x is either due to the formation of Mn-rich spa- 
tially ordered regions, or to n-p co-doping effects from 
the n-type oxygen impurities, or a combination of both. 
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